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' . THEORETICAL STUDY OF AIR FORCES ON AN OSCILLATING OR 
STEADY THIN WING IN A SUPERSONIC MAIN STREAM 
By I., E. Garrick and S. I. Ruhinow 


•SUMMARY 


A theoretical study. Based on-the linearized equations of motion 
for small disturbances, ie made. of the air forces on wings of general 
plan forms moving forward at a constant supersonic speed. The 
Boundary problem is set up for Both the harmonically oscillating and 
the steady conditions . Two types of. Boundary . conditions are 
distinquished, which are designated "purely supersonic" and "mixed 
supersonic" . The purely supersonic case involves independence of 
action. of the upper and lover surfaces of the airfoil and the present 
analysis, is mainly concerned with this case. A discussion is first 
given of the fundamental or elementary solution corresponding to a 
moving source. The solutions for the velocity potential are then 
synthesized By means of integration of the. fundamental solution for 
the moving source. The method is illustrated By applications to a 
number of examples for Both the steady and the oscillating cases and 
for various plan forms, including swept wings and rectangular and 
.triangular plan forms. The special results of a number of authors 
are shown to Be included in the analysis. 


INTRODUCTION 

This paper constitutes a theoretical study of the aerodynamic 
forces on an oscillating or steady wing of finite span moving forward 
af a uniform supersonic speed. . The treatment is Based on the 
linearized theory obtained By considering only small disturbances in^ 
an ideal fluid. The wing is therefore considered to Be a nearly flat 
thin surface at a small angle of attack and the flow is considered 
nonviscous and free of strong shocks : .- The theory in this case is 
equivalent to finding .certain solutions of the. wave equation in three 
dimensions with respect to. a moving coordinate system. 

For the case of steady motion there exist a. number of interesting 
solutions and methods . Among these may be mentioned the von Karman and 
Moore, linearized treatment of slender Bodies of revolution (reference 1 ) 
the Prandtl acceleration— potential method; employed - by Schlicting 
(references 2 and 3) > the Busemann.method of "linearized conical flows" 



2 


NACA TN No. 1383 


(reference 4), studies of Jones, Puckett, Stewart, Brown, and 
Gurevich (references 5 to 9); and a method of von Harman employing 
Fourier, integral solutions of the two-dimensional wave equation and 
described by him as "acoustic oscillator method" (Wright Brothers 
Memorial lecture, Dec. 17, 19^). 

The corresponding unsteady or nonstationary problem for two- 
dimensional flow (infinite aspect ratio) may be considered to be 
solved. In this connection there may be mentioned the work of Possio, 
von Borbely, Temple and Jahn, and the present authors (references 10 
to 13). Of 'interest also are . two wartime German papers by Schwarz 
and Honl (references lk and 15) • The corresponding steady plane 
case to which the nonstationary problem may be reduced is that 
treated by Ackeret. ' 

Eesults for the nonstationary or oscillating case are of great 
interest in the investigation of aircraft instability. The two- 
dimensional results have been applied to a study of flutter at 
supersonic speeds, in references 12 and. 13- Of more direct interest 
for this application are the three-dimensional results, especially 
for wings of swept plan form. 

The method used in the. present study is to build up solutions of 
the equation satisfied by the velocity potential by superposition of 
the fundamental wave -potential solution for a spherical source. These 
solutions are also made to satisfy certain required boundary conditions 
on the airfoil surface . In the two-dimensional supersonic nonstationary 
case, which appears herein as a '•special limiting case, ft can be 
proved that the procedure leads to a solution that is the unique 
solution of the given . boundary problem. (For the problem of subsonic' 
flow past a thin wing reference may. be made to the general treatment 
and method of Kiissner (reference l6) which also involves solutions of 
the wave equation.) : ' 

Some qualitative features of the nature of the boundary problem 
may be mentioned here,. Further remarks may be found in reference 17 
and in Von Harman *s Wright Brothers Memoria.1 lecture. In the case of 
subsonic flow past an airfoil' the whole field is influenced by the 
body. The concept of circulation has proved to be very useful and the 
Kutta condition has been .used to specify the circulation by requiring 
smooth flow leaving the trailing edge- Thus, a deflected aileron in 
subsonic flow influences the flow pattern over the whole wing even 
more importantly than over the aileron itself. ; 

In the case of supersonic flow the influence of the body is 
limited to only certain parts of the field of flow and generally the 
wake does not influence the upstream flow. The boundary problem for 
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a three -dimensional surface moving at a supersonic speed can "be 
classified into two types, referred to herein as "purely supersonic" 
and 'inixed supersonic The definition of these terms Is given in 
the analysis according to the parts of the field influenced "by the 
airfoil, the purely supersonic case involving independence of action 
of the top and "bottom surfaces and no reflecting surfaces in the field. 
Thus, in the purely supersonic case, a deflection of the aileron 
would produce only a local effect at the aileron; in the mixed 
supersonic case, it may have a decided influence on the part of the 
wing adjacent to the aileron or on other parts of the wing. For a 
given wing "both types of problems may be involved. 

The treatment used for the purely supersonic cases, involving 
source and sink distributions to account for the action of the body, 
is believed to be exact within the framework of the linearized theory. 
The upper and lower, surfaces of the airfoil are regarded as acting 
independently, each surface being ''unaware" of the presence of the 
other. The treatment is thus analogous to that of sound in a moving 
medium generated by the motion. of pistons imbedded in an infinite 
plane. This flow picture is obviously incomplete in -the mixed case 
and more complicated di8tributions (doublets) are also required. 
However, for some purposes, the simple? 1 treatment may still be used 
in conjunction with appropriate correction factors. Also, for . 
steady flow past a symmetrical airfoil at zero lift, the simpler 
treatment, oo.n be employed for study of the wave drag. 

The object of the present paper is to develop the expression for 
the velocity potential in the purely supersonic case, based on the 
elementary solution for the sound scarce moving unifonrtly at a 
supersonic speed, and to indicate it 3 application by a number of 
special examples . - . 


ANALYSIS 

Wave Equation and Source Solutions 


. In the. linearized theory based on small disturbances; the equation 
satisfied by the velocity potential for the propagation of sound waves 
of small amplitude is the wave equation 



4 


NACA TN No. I383 

Symbols are .defined in appendix A. . The fluid medium is considered 
at rest at infinity. 

In the treatment of linear, partial differential equations the 
so-called elementary or fundamental solution is of great importance 
since general solutions cari he "built up by distributions of 
elementary solutions . From a physicd]. point of view the elementary 
solution may correspond to a source. A discussion of the nature of 
elementary solutions for hyperbolic differential equations of a 
general type has' been given by Hadamard (reference 18), who makes 
the cardinal statement that "every result of the theory can be and 
has to be deduced from the, consideration of the elementary solution 
only". ... ... 

A fundamental solution of equation (l) from which general ’ 
solutions may be formed is that of a source of sound fixed in the 
medium ' 



c- 


where ", 

" r ’;= !/(? - J'f >;(y' - V) 5 + (z’ - 

In equation (?) the fixed source is located at the point (f/, tj', t,'), 
the strength of the source is A( q £')f(t'), and the minus . 

sign indicates 'that the spherical waves are diverging from the . 
center of the disturbance. 

Another closely related solution of equation ( 1 )' is that of a 
fixed point source for which the spherical waves are converging onto 
the source 

A-;; 4' *.t) ; 

The wave potential in equation (2) is often designated "retarded" 
and that in equation (3) "advanced". 

It is intended to consider thin lifting surfaces of small curva- 
ture which are moving forward at a constant supersonic velocity v 
and which may be performing small oscillations normal to the direction 
of v. The direction of v will be that of the negative x-axis and 
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the surface will he replaced by a distribution of moving sources 
in the plane x, z (fig. 1). 

Consider a source moving in the negative x direction with 
uniform velocity v and a rectangular coordinate system attached 
to the moving source. If the new coordinates are designated 
hy x, y, z, t, where x = x' + vt', y = y', z = z', t = t', the 
equation satisfied hy the potential is 



or 

i_ 2v Af \ h 2 (f> c!y 

c 2 St 2 c 2 bxbt Vc 2 ' ) bx 2 by 2 bz 2 

This equation is satisfied hy the potential of sources of sound in 
motion through the medium with uniform velocity v in 'the negative 
x direction i It is also the equation satisfied hy the disturbance 
velocity potential for a fixed body creating small perturbations 
from an oncoming main stream of velocity v in the x direction. 

A brief derivation from hydrodynamical principles is given in 
appendix B. 

It is known from the classical study of the wave equation 
(reference 16) and can he verified hy direct substitution that a 
solution of equation (1) is transformed to a solution of equation (4) 
by means of the following substitutions, corresponding to a combination 
of the Lorentz transformation and a Galilean transformation: 



where M, the Mach number of. the main flow, is v/c. 
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For the purpose of studying the supersonic case (M >1), it 
is more convenient to employ modifications of transformation (5) 
obtained by multiplying the right hand Bide by the constant 

1/ l/l - M 2 , or . 


x‘ = 


y = 


z' = 


1 - M 2 

y 

/l - M 2 
z 

l/l - M 2 ' 


t* = t + 


xM 


c(l - M 2 ) 


(5*) 


The particular solution of equation (4) that corresponds to a 
moving source will be seen in the following discussion to be analogous 
to a solution of equation (l) given by the sum of potentials in 
equations (2) and (3), namely to 


*•’$(*' - 7 ) * 4' * r)l 


( 6 ) 


The desired solution of equation (4) corresponding to equation (6) 
is -obtained' with the aid of the substitutions (5a) as ■ . ' 


<k - 


A 


•L - - - £\v r(t ■ 

\ . - :c M 2 1 °/ \ c m - 1 <V 


(7) 


where 


r = 


M 2 - 1 


1 /(X -s) 2 - <:r - •-) [(> - 1 ) 1 ' + (Z -t V j (8) 


(The term ' in equations •(53) causes no difficulty since 

only the squares of the space coordinates are needed.) 
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This solution for may be expressed in the form 


where 


*0 ■ r[ f (* • T a) + f (* • T l)] 
M x ” - r 

M X ’S _ r_ 

1 ' c M 2 - 1 c 


(7a) 


and where r is defined as in equation (8). The constant A (c., q, £ ) 
could of course have been included in the functional symbol f but has 
been separated for convenience. It may be considered to represent the 
space variation of the source strength as distinguished from the time 
variation of strength. For a moving source of constant strength the’ 
time function may be considered equal to unity and the potential 
expressed as (referenced): 



It will- be recognized that the solution, equation (7a), is valid in 
a conical region, the so-called "Mach cone", opening aft- of the 
moving source . Outside of this conical region, defined by the 
equation r = 0, the flew is undisturbed- 

The result expressed by equation (7) may be considered physically 
from two points of view. In one, as considered by Prandtl (reference 2 ), 
a source of variable strength moving along a certain path iB replaced 
by a continuous succession of fixed source pulses distributed along 
its path acting consecutively one after the other. Each pulse, . 
considered fixed in an absolute coordinate system, emits a spherical 
wave traveling at sound speed, and the coordinates of the center of 
the spherical surface are £ + vt, q, £. The radius vector R of a 
point x, y, z with respect to this center is 

R =V[x - (£ + vt)] 2 + (y - q) 2 + (z - t ) 2 
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The time at which the spherical wave passes the point x, y, z is 



c 


Eliminating R between the preceding two relations results in 


c 2 t 2 - (x - vt) 2 - (y - t|) 2 - (z - £) 2 = 0 


The roots of this quadratic equation in t . are precisely the 
quantities and t 2 defined in equation (7a) . That is, the 

field point (x, y, z) is influenced at time t by two waves which 
originated at times and T-| earlier • It is of interest to 

observe that, in the supersonic flow, both roots are real and positive 
and have physical significance; whereas, in the subsonic flow, only 
one root is positive and of physical significance. In the supersonic 
case the field of influence of a source is the particular Mach cone, 
with vertex at the source and, through each point in this region at 
instant t, there pass two spherical surfaces representing the waves 
originating at times and t 2 earlier (fig. 2). 

From the other point of view of the result (equation (7)), a 
single diverging spherical wave -pulse is considered. Let this wave 
originate at the point (£, 11, t ) at a time T (fig. 3) and consider 
its effect at a point (x, y, z) (within the Mach cone whose vertex is 
at (£, tj, £)) moving with a velocity greater than that of sound. 

Clearly at a later time (T .+ t^) the moving point penetrates the wave 

front and at a still later time . (T + Tg) it emerges from the wave front. 

The potential at . x, y, z changes only on entering and on leaving 
the wave front and the two terms in equation (7) correspond to these 
two effects . The factor ■ 2 appearing in the potential for a constant 
source moving at a supersonic speed also has its origin in this 
physical fact, in contrast to that for a source moving at a subsonic 
speed, where the field point penetrates the wave front but never 
emerges and where the corresponding fax tor is unity. The two- 
dimensional supersonic case involves cylindrical waves and the potential 
of the point x, y is continuously changing from the time the point 
enters to the time it emerges from the wave (reference 13 ). Observe 
the interesting geometrical property of r (equation (8)), namely 2r is 
the difference of the radii of the spherical wave at time t 2 and at 

t ime T]_, that is, r = ( T 2 - T 1 ) . (Observe also that the potential 
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which formally appears in equation (6) as the sum of potentials, 
half -advanced and half -retarded, transforms in the moving 
coordinates to a sum of retarded potentials in which the original 
retarded part is associated with the diverging spherical concave 
wave from which the point is emerging and the original advanced 
part- is associated with a diverging convex wave into which the 
point is penetrating.) Recent papers of interest in connection with 
moving acoustical, sources are references 19 and 20. 

Surface Distribution of Sources 

Sources and sinks of the type; 0 O will now be distributed to 

represent the upper and lower surfaces of a thin airfoil. The 
procedure to be followed is that used in the two-dimensional case 
.(references 10 to 13) where the upper and lower surfaces are 
considered separately. Also the total effect may be separated into 
an effect of the mean camber surface and an additive effect due to' 
thickness alone. In most of the applications, unless stated to the 
contrary, the mean camber surface is- considered. 

Let a continuous distribution of sources be given over the mean- 
camber surface. The airfoil is considered so thin and flat that the 
source distribution may be treated in the x, z plane (fig. l). The 
airfoil" surface may be considered moving at a constant speed v in 
the negative x direction (or fixed in a stream moving in the 
x direction). The effect at a point x, y, z at time t of a 
distribution of sources of position magnitude A(g, 0 , £ ) i 3 given by 
an appropriate integration over a region of the , f plane of the 
form ' 1 ' ' . “ ' 


f(x, y, i, t) = J I <$ o d£ d£, ( 9 ) 

where represents the function given in equation (7) with rj = 0 . 

The total effect at the point (x, y, z) is the sum of the effects 
of all disturbances having their origin within the Mach cone with 
vertex at x, y, z and opening in'- the upstream direction. This 
conical region need not extend into the undisturbed part of the flow, 
that is, it need not extend beyond the most forward surface' envelope of 
the Mach cones of influence of the body. There are essentially two 
types of boundary conditions that need to be distinguished, designated 
by the terms "purely supersonic" and 'mixed supersonic". A point of 
the boundary belongs to a purely supersonic case if the upstream facing 


10 


NACA TN No. I383 


Mach cone contains, in the part of the x, z plane not considered 
occupied by the tody, no disturbed fluid having a component normal 
to ' this surf ace . Otherwise the point belongs to the mixed supersonic 
case. A sufficient (hut not necessary) criterion for. the purely 
supersonic case is that the component of the main stream normal to 
any edge or contour of the plan form in the x, z plane (contained 
within the upstream facing Mach cone of the given point) shall he . 
supersonic. There is no downwash ahead of the body, no holes are 
in the body, no spilling of fluid occurs around edges, and no 
reflecting surfaces are in the flow field. In this case the upper 
and lower surfaces of the airfoil, ere considered to act independently 
of each other j a disturbance created on one side does not affect the 
opposite side.. The flow can be considered to arise from the 
appropriate movement of small' pistons acting at the regulating or 
generating surface. This, condition is in contrast to that of the 
mixed supersonic case, for wjiich the effect of the disturbance spills 
over the edges or sides, and a disturbed fluid region (downvash) 
may exist. ahead of the body. Thus, points of a triangular surface, 
moving vertex foremost and completely, outside of the Mach cone 
associated with the vertex, belong to the purely supersonic case. 

If the triangular surface, is inside the Mach, cone associated with 
the vertex, the points belong to the mixed supersonic case. Of course, 
for a given surface, both cases may be involved. A few examples are 
shown- in figure 4 . _ 0 

In the purely supersonic cane the circulation concept plays no 
particular role and the drag associated with lift or thickness may 
properly be denoted as wave drag. In the mixed case the flow retains 
subsonic features and the drag associated with the lift is sometimes 
denoted as induced drag. 

Although the treatment given for the purely supersonic case is 
believed exact within the limitations of the linearized theory, an 
exact treatment of the mixed supersonic case is not available • These 
problems involve greater difficulties in the boundary conditions, for 
the flow to a certain extent acquires features of a subsonic flow in 
that the fluid field "senses" the approach of the body’. Thus, in 
certain cases conditions at the leading edge, at the trailing edge, 
and in the wake must be specially taken into account'.- However, for 
some purposes and in certain prob.lems it may be useful to treat the 
mixed supersonic case in the same manner as the purely supersonic case 
and to introduce appropriate correction factors • 

The region of integration in equation (9) is the part of the 
body (in the £ plane) cut but by the upstream opening Mach cone, 
with vertex at~(x,' y, z). This region in general depends on the plan 
form of the body as well as on.(x, y, z). With the understanding 
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that the leading point of the. "body is at = 0, the integration 
may he written 



where 



•The limits of integration £^ and £ in equation ( 10 ) may he 

recognized as the distances from the £ -axis to the near and far 
sides, respectively, of the hyperbola defined by the intersection 
of the cone r = 0 and the plane q = 0 . Thus , from equation ( 8) , 
with t) = 0, and £0 are recognized as the roots of the 

equation 

iiFWk 7 ?) - 0 . < u > 

The limit ^ in equation ( 10 ) represents the § coordinate of the 
vertex of the hyperbola and is defined by the condition ' ’ £ = £ 

that is by £^ = 0 . The point (^, z) is the farthermost downstream 
point which can affect the point (x, y, z). 
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Boundary Condition 

The strength of the distribution of singularities in equation (10 ) 
will now "be determined "by the "boundary condition of tangential flow 
along the airfoil surface. The "boundary condition may be expressed as 

(~\ - v(x, z, t) 

, + ' ( 12 ) 

8x $t 

where the airfoil shape is defined by y = ,g(x, z, t) and where the 
two terms represent the normal Velocity induced by the airfoil shape 
and by its own proper motion. It is shown in appendix C and can s.lso 
be made clear by physical reasoning that as' y approaches zero from 
the positive side (y --■ + 0) 

M- = -2n(ir - 1) A(x, 0, z) f(t) 

Sy 


or, briefly j 

A(x, z) f(t) = - i— w(x, z, t) (13) 

2it(rr - 1 ) 

As y approaches zero in the negative half plane., ah equal and 
opposite result is obtained. Equal source distributions on the upper 
and lower surfaces therefore result in a discontinuous vertical 
velocity distribution near the plane y = 0 and may be used to 
represent symmetrical, thickness distributions .• "The source distribution 
representing a thin body with arbitrary thickness distribution is in 
general unequal on the two surfaces. The effect of thickness is 
discussed in a separate section. A representation of the mean-camber 
surface alone may be obtained by placing equal and opposite sources on 
the under surface in proximity to the sources on the upper surface. 

The potential (ft is to be understood in the subsequent analysis to 
be prefixed by a 1 sign, plus for the upper surface and minus for the 
lower surface . The vertical velocity will in general be measured 
positive upward. 




where, for y = 0, (see equations (7a), (10), and (11)) 

* *1 ” * ~ 1 - 1 ■ " * ■ ■ 


. cp 2 . . cp 
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and where it Is understood that W(£, £) = 0 at any point off the 
"body or .where the .integrand is not real. 

Equation (15) may W put into a - simpler 
of a hew variable 0 instead of £, which is obtained from the 
relation (see appendix C) 


H = (^2 " £l) cos 0 + + ti 


or 


t = ? cos 0 + z 

b bQ 


The surface potential (equation (15)) may then be written as 


nx nn 


Jf(x, z, t) = - 


r 


*** Jo 


/ 0 


w[§, £(0)|jw(t - + w(t - T^] d0 d§ 


(15a) 


where 


* ■§ 


T 1 = 2 

7 cp- 


- (M - sin 0) 


:: ' X - £ 

t = - — — (M + sin 0) 
2 cp 2 

-1 ^ ‘ Z 

0 = cos P . 

' * - £ 


I . 

. Equation (15) represents the central result of the analysis and 
within Hie limitations already discussed may be applied to wings of 
any plan form in steady motion or performing small’ oscillations . In 
the stationary or steady case, (j> does not depend on time and the 
function w(t) is to be replaced by unity- Then, in equation (15), 
wYt - t \ + w/t - is to be replaced by 2. 


/ 
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Pressure Relations 

For tiie sake of reference, relations for the pressure and the 
lift and drag forces are given here. The disturbance pressure 
(local static pressure minus the pressure in the undisturbed stream) 
may he written as 





(16) 


The pressure difference (positive if acting' downward) at any point 
(x, z) may he expressed as 


% • Pl 


where the subscripts U and L refer to the upper and lower surfaces. 
For the mean camber surface Pj. = -p^ and 



i 


The total. forces on the airfoil in the .. y and x directions are 
given by ' ’■■■■' 


Y =, Lift =. 
X = Drag as 


//« 

si" 


da 


dy dz 



where the integration is to be taken over the 
Expressed as integrations over the plan form 

Y = JJ (*L " %j dz 


r n f /dy\ . _ /dy\ 


complete airfoil surface. 


( 18 ) 
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It Is often convenient to separate the . slope terms as follows 



where a is the conventional direction of the main stream with 
respect to a reference chord, and and . <j^ are the local slopes 

of the airfoil surfaces measured with respect to the reference chord 

and positive in the same sense as a. 

! 

APPLICATIONS 


Wing of Infinite Span and Zero Sweep 

For the first application of equation (15) the results for Doth 
the oscillating and steady two-dimensional case will "be derived. For 
the harmonically oscillating wing having identical motion in every 
chordwise section, the vertical velocity can he written in the 
complex form .... 


w(x, t) = W(x) e iayb 


Then 


w 


> . to 4 imt / — icoTp — InyT'A 

(t - t\ + w^t - = e^° (e 2 + e L J 


( t 2 +t 1 

. 2 


T 2 “ 


COS CD — 


T ~) 


Equation (15 a ) becomes 


(x, t) = - 


icot fix -ia r^ /** / x " S cd. 

e 1 — . cfi d I cos gin e j d0 ^ 


e icot fix 

Jo 


W(£) e 


3 


’) 


(19) 


where 3 


, f¥ 


1. 
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The integration with respect to 0 may "be readily performed 
with the- aid of the relation 


Finally 



cos (\ sin 0) d0 = J D (X) 


e ipt Ox 

Hx, t) = g- W(£) I(£, x) d £ 

do 


( 20 ) 



This result for the velocity potential is identical with equation (ll) 
of reference 13 and is used therein as a las is for calculation of the 
nonstationary two-dimensional case. 

In the steady case, p' = 0 and !(£, x) = 1. The expression for 
the velocity potential is 


0 (x) 



x 

W (|,) d£ 


where 


W(£) = v — . 

dx 


This formula or the pressure relation 


P = 




dy 

dx 


( 22 ) 


applied to loth the upper and lower surfaces of the airfoil leads to 
all the results of the Ackeret theory. 


Wing of inf ini te' Span with Angle of Sweep 

• Consider an infinite wing with angle of sweep A (fig. 5)-/ and ■ 
assume that all sections in the flight direction are identical in 
shape and that the wing is undergoing harmonic motion ^ In general 
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the vertical velocity w can he written in the complex form 


w(x, z, t) = W(x, z) e J 


If each section normal to the leading edge is performing the same 
motion the form of W(x, z) is W(x - z tan A). If the wing is 
assumed to perform pure vertical motion alone, then W(x, z) is 
a constant. If the wing is assumed to rotate about an axis 
x = constant, then W(x, z) is of the form W(x). 


The potential is of the form (fig. 5 ) 


0(x, z, t) « 



F de d£ 


4 P e i 



F d 0 d£ (23) 


where 


1 


Ffe, e, £(»).) Je“V V 




and where 


83- 


ZB 


1 - (3 cot A 


£4 


X + -Ufj 


1 + 3 cot A 


e, - 


-1 

- cos 


^CQt 


A 


The values of the limits £. and are found hy solving for £, 

in the relations = Z, cot A and £ = £ cot A which represent 

the intersections of the Mach lines through x with the leading edge 
The limit 0 = 0 corresponds to £ '= £ cot A, the leading -edge line 
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When ¥(.*, z) is a constant or a function of £ only, the 
velocity potential can he expressed as - 


(/>(x, z, t) 


e tot 

P 



w(£) i(£. 


x) ag. - 



W(S) I.^, x, 


*). d£ 

(24) 


where I (c . x) is as defined previously and 


n ( J> x, z) , i. too '-'e? J o cos a sin 


ae (25) 


Observe that the integral involved in equation (25) for I 1 reduces 

to the Bessel function of zero order when Qj_ = n as in equation (l8). 

Tills interesting integral may . therefore he called an "incomplete" 
Bessel function of zero order. Systematic investigation of its 
properties would appear to he desirable . 


For the infinite swept wing in the steady case the frequency o> 
may he made equal to zero in equation (23). Consider as a simple 
example the case of a thin wing at a small constant angle of attack a, 


dy 

that is, — - -a. let the angle of sweep he less than the complement 

U-A 

of the Mach angle, that is, p cot A >1 (otherwise the case involves 
the mixed -supersonic flow conditions). From equation (24) with a* = 0, 

. / x f' n 

1(2,, x) -- 1, and !]_(£,, x, z) = — — ; , 


r /('h = 


va 


d£ - 





cos 


-l/Scot A - z \ 

(-~7 j as 


• x cot A - z 

.== va 


\/p 2 cot 2 A - 1 


(26) 



20 


MCA TN m .. I383 


The local pressure difference is given "by 

O 

2pv a p cot A 
0 |/| p 2 cot 2 A - 1 



This equation reduces, for 


A = 
*o 


0, to the A ckeret result 

- g P v2(X 
P 


( 27 ) 


Let the index n refer to quantities measured normal to the leading 
edge . Then 

oc^ = a sec A 
. v n . = v cos A 

= ~ = M cos A 

and 2 

_ g <”n °n 

p " rtt 2 ■ 1 

a result similar in form to the expression for p D and already 
stated by Busemam (reference 17) in 1935- (See also reference 6.) 

The harmonically oscillating case with W(x, z) assumed to be 
of the form W(x - z tan A) leads in a similar manner to a result 
analogous to equation (20). 


Rectangular Wing of Finite Span (Zero Sweep) 

Consider a harmonically oscillating rectangular wing of finite 
span as in figure 6. Region I is described as purely supersonic and 
region II as mixed supersonic. The higher the aspect ratio and the 
stream Mach number, the relatively smaller the region II becomes . 
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The potential for region I for identical motion of each chord - 
wise section is exactly that given for the infinite wing in 
equation (20). However, more general types of motion involving 
spanwise variation may also "be treated-. For example, let . the wing 
perform harmonic oscillations in vertical tending and in torsion t 
about a spanwise axis x = x Q In certain prescribed spanwise 

modes . Then, with a and h used -to describe angle of attack 
and vertical position (fig- 7) 



where the various tf's are associated with the corresponding 
variables in equation (29). 
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With the use of equations (15a) , these potentials maybe . 
"expressed as 



where 


1 



- l) 

c(^ - 1) 


and, expressed as a function of 9 3 


■ - OqC 2 + £ 0 cos 0) 

. h i(0 = Mz +> t Q cos e) 

I 

If the modal functions in equations ( 31 ) are a, = h^ = 1 , the 

potential corresponds to that given by equation (20) for the two- 
dimensional case (see also equation ( 14 ) of reference 13 ). It is 
of interest to consider modal functions for and h-j of the 

type 1 (t/s) n where s is the semispan. For modal functions of 
this form the- typical integral involved in equations . (31) may be. 
expressed as 
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With the substitution of - 0 for 0, F n may be written as 


,TT 


n ~ J o ( z + t 0 ein s ) cos ^ cos 0'j dff 
( z - sin 0) cos cos 0^ d0 


The further reduction of F n is cade with tlie aid of the following 
relation (reference 21) 

J(x , n ■ ' * 

jk'*' 7——; — "x I cos (^ cos 0 ) sin 0 d0 




?? e , C tf e n " 0 corresponds to constant modes and yields 
for the potential the result already given by equation ( 20 ). The 

case n - 1 corresponds to linear modes and the function F x becomes 


F, — tczJ . 

1 o\m 


^ U ^ ll2ed in the equation for the potential yields £ 

ih 8 i t 18 rm Ule two -dimensional -case result multiplied by the 
factor z/s. The case- n = 2 corresponds to parabolic modes and 


the function Fg becomes ? 


.-gj 


When F 2 is used in equations ( 3 l) } the J Q term yields an integral 
of the type given by equation ( 2 0 ). With, the use of the relation 
l lA; " _J o tae J i term also yields an integral of the same 

form.. This type of reduction to the form of equation ( 20 ) my be 
made in general for any integral index n by means of the recurrence 
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formulas for Bessel functions^ and thus use may he made of the 
numerical procedures used for equation (3.0)* (See reference 13*) 

It may he of interest , to treat the potential for the mixed 
supersonic region II (fig. 6) as. though it were part of a purely 
supersonic region. The equations corresponding to equations (31) 
are „ 



f e -H 

Jo 


°i (0 


cos 



sin 0] d0 


•) 






Oq_(£) cos j^-sin 0^ . d0 d£ 


(33) 


and similar equations for: ^ and (f^- The limit ^ is found as 
the' value of £ for which \ = s or 


gg « X - (8 “ Z)P 

The limit 0=0 corresponds to £ = £ . and the limit. 0 ='0 2 
corresponds to. £ = s , or, from equation (15a), 



The last term in equation (33) leads to integrals of the "incomplete". 
Bessel function type as mentioned for the case of the infinite wing 
with angle of sweep. 

The foregoing results for the oscillating rectangular wing will 
now he specialized to the steady case (o> = 0, q = 0, a^(£). = 1, 

a^(t) = a, the constant angle of attack) . Then, from equations (31), 

the velocity potential for region I # is 



( 34 ) 
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For. region; II, from equation ( 33 ).; , : ; . i : ' • 



The actual integration in equation (35) may he easily perforated 
"but is hot required for the purpose of obtaining the local 'pressure • 

The local pressure difference is directly obtained for regions: I 
and II from equations (34) and (35), as 

V 

2 • -• \ \ 

'■ 2pV a " •• - : 



tip 


It my he observed that p is constant along rays from the 

g • 2 

= Constant. Along the ray corresponding to the Mach line 

from the tip - — ■ — 3=1 and p takes on the constant value p . 

X • .. I 

Along the ray corresponding to the tip 2 •= n, half of this value, 
is obtained* This edge condition is physically incorrect since the 
assumption of the independence of the two surfaces of the airfoil is 
not correct near the tip... . . 


This particular problem has been treated by. Busemann (reference 4) 
by his method of conical or perspective symmetry. The condition 
along the ray corresponding to the tip is p = 0 and Busemann' s 
result for region II is 


2pv 2 a 1 


•II 


it 


cos 




2(3 




The total lift over region II is one -half of that of an equal area 
of region I. A comparison of this result and equations (36) is 
shown in figure 8. This comparison gives an indication of the 
errors involved in the assumption of independence of the two surfaces 
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near the rectangular tip and conversely it gives an indication of the 
appropriate correction factors required to allow for the tip effect . 
It appears that equations (36) overestimate the lift over all of 

2 /* 

region II "by a factor 1 - — or "by approximately 3 ° percent. 

Thickness Distribution 

It has already been remarked that the treatment employed in 
the analysis mainly for the mean-camber surface can also be applied 
to obtain the effect of thickness.' In equation ,(15) the vertical 
velocity ¥(£,, £ ) may be specified for both the upper and the lower 
surface . 


As an example, consider a plan form such as that shown in 
figure 9 in steady supersonic flow. Let the airfoil section shape, 
for convenience chosen symmetrical and independent of span, be 
defined in the center section by y = g(x) (and in any other section 
by y = g(x - z tan A)). Then, for the upper surface, 

W(x, z) = va + vg 1 

and, for the lower surface, 

W(x, z) = va - vg' 


where g ' is. the derivative of g with respect to- its argument . 

The velocity potentials in the various regions in figure 9 are 
of the form 






f ae as,. 




I»§> 




oe r 


09, 


F d0 d£ 


+ f 


/I0, 


no nit 

k 3 k: 


*> 

‘ fir 


F d 0 d£ 


>( 37 ) 


F d 0 d c, + 


JO (/0, 


F d 9 dc, 


J 
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The limits in the foregoing integrals are as follows: the 

limits 9 = 0 ^ and 0 » correspond respectively to the leading 

edge lines £ =£cotA and £ = cot A; 0' = 0 and 9 - n 

correspond respectively to the Mach lines £ = £ 2 8114 £ = 

and are obtained, respectively, from the relations 

= £ cot A and = -J, cot A; 0 = 0 2 corresponds to £ = s; 
is obtained from the relation £ 2 = a; 5 ,^ is obtained from 
£g = £, cot A; and is the value of £, for the leading edge 

of the tip. Then 


0t = 



-ll 

cos ! 

cot A - z 

’1 - 

^ at - g 


-1 

cos 

/- c, cot A ■ 

\ X - § 


-1 

cos 

s - z 

fi ■ 

'2 = 

X 


") 


6 , = 


6 « « 


f _ X • Zp _ ¥ X + Zp^ ^ X “ Zp 

h " 1 - A ' b i'r " 1 + A ' ^6 " 1 > A 


=.x - (s - z)p ; = a tan A ; A = p cot A 


If, for example, the distribution function F is a constant K 



x cot A - z 
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The corresponding local pressures are 


piii 


pvK 

A 

3 

\A a - 1 

pvK 

A 

3 

ro 

1 

H 

pvK 

A 

3 

\Jk 2 -1 

pvK 

A 

3 

^ - i 


u 1 -11 + ab 

1 _ ± COS 

V * A + B 


^ ( J 1 + AC ' r . -1 1 - AC 
=■ ( COS — + cos : 


A + C 


A — C 


2 _i k c - 1 
- tan \ /— — - 

1 - C 2 


where 


A = 3 cot A ' 

B = - z ) p 
X — s tan A 

C = - 3 ' 
x 

The constant K may he interpreted as va associated with 
constant angle of attack. In this case, region II is to he regarded 
as a mixed supersonic region and the result given is not the 
appropriate solution for this region. If the constant K is 
interpreted as g* = Constant, the results are applicable to a thin 
symmetrical wedge of half vertex angle K and may he employed to 
yield the wave drag according to the linearized treatment. 

Jones (reference 5) treats symmetrical airfoils of various plan 
forms at zero lift hy use of pressure potential. The use of velocity 
potential leads to the same results as given in reference 6. Thus, 
equations (13) and (14) of reference 5 for a wedge correspond to the 
preceding results. The velocity potential in general is more useful 
to treat pressure distributions for a given body whereas the pressure 
potential may be more readily adapted to treat airfoil shapes and 
plan forms associated with desired types of distributions of pressure 
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Triangular Plan Form 

The triangular wing (fig. 10 ) extending across the Mach lines 
from the vertex may serve as a f inal example. ^ For the steady case 
of a vanishingly thin surface at angle of attack a. the velocity 
potentials and pressure relations for regions I and III are 
equivalent to those just discussed in the preceding section. The 
lift AL on a strip Ax of the triangle located at abscissa x 
from the vertex is given by 

px cot A 

AL = Ax / Ap dz 

J-x cot A 



A ,1 

- 1 ♦ 1 . | /V - *) 


L p 2 • 

A + 1 p 2 \ ' A + l/_ 


where the two terms correspond to the integrations over regions I 
and III, respectively. Then, 

AL = — ■ Ax Ax ■ 

M 2 - 1 


The area of the strip is 2x Ax cot A and hence the lift coefficien 
is independent of x and equal to 


AL 



4a 

P 


Gurevich (reference 9 ) treats this case and his relations can be 
shown to be equivalent to the foregoing ones. The pressure 
distribution is illustrated in figure 10, where p Q , the reference 

pressure, is 2pv 2 a/p. Observe that the pres sire area above the 
unit ordinate cancels the area of pressure deficiency below the unit 
ordinate. Also shorn in figure 10 is the distribution of pressure 
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» 

as the half vertex angle of the triangle approaches the Mach 
angle • 1 - 

The triangular wing inside the Mach cone from the vertex 
requires a more elaborate treatment (references 7, 8, and 9). 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va-, June 4, 1947 
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APPENDIX A 
SYMBOLS 


4 . 
to 

X', y', a* 
x, y, z 

S', V, £’ 


l* ^ t 


t, T, f 

V 

c 

M 

r 


T l' t 2 
A 

p = <?" - 1 

g 

t 0 , t 2 . Si 


9 


disturbance -velocity potential 

potential function defined in equation (7) 

rectangular coordinates for fixed system 

rectangular coordinates attached to source moving 
in negative x direction; also represents field 
point being influenced 

rectangular coordinates used to represent space 
coordinates in fixed system 

rectangular coordinates used to represent space 
location at source distribution A( q, £ ) 

time 

velocity of main stream 

velocity of sound 

Mach number (v/c) 

distance defined by equation (8) 

time function defined in equation (7a) 

in applications, (3 cot A 

function defining airfoil surface (y = g(x, z, t)) 
limits defined in equation (10 ) 

variable used instead of t defined by relation 
preceding equation (15 a ) 


P 


pressure 
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Po 

P 

a 

a 

O) 

<1 

v(x, z, t) 

A 

h 

h 


reference pressure 
density 

angle of attack 

time derivative of a 

angular frequency 

relation defined in equation (33.) 

vertical velocity factored in equation (lk) as a 
space function W(x, z) and time function w(t) 

angle of sweep 

vertical displacement . 

time derivative of h 
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APPENDIX B 

DIFFERENTIAL EQUATION FOB TEE VELOCITY POTENTIAL 


A derivation of equation (L) is given "briefly here. The 
condition for irrotational flow is 

curl "v = 0 (Bl) 

and this relation implies that a scalar velocity potential ft 
exists, such that 


. v = grad (j 

The general equation of continuity 


— — + div pv = 0 
8t 


(B2) 


may he written as 


il ?E + v 2 0 = 0 

p Dt r 


(B3) 


where differentiation f05.l01ri.ng’ the particle is denoted "by 

d a ._ 

grad) 


and v = div grad is the Laplacian operator. 

From Euler 's equations, or from the general Bernoulli relation, 


v 2 

T — + + 

dt 2 


^-0 


(bL) 


where a space constant function of time has been included in and 
where it has been assumed that p is a function of p only . 

With the use of equation (BL) and the acoustic relation. 


2 _ £p 
= dp 
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where c is the local variable speed of sound, it follows that 

c 2 , 

= -—grad p 


l;Sp 

P St 

£ Sp 

p St 

With the aid of these two relations the first term in equation (B3) 
becomes 

iSiL. . grad 

P Dt c 2 ^ t 2 St 


For sms.ll perturbations from the main stream of velocity v in the 
x direction, c may be considered, equal to the constant speed of 
sound in the undisturbed medium and, in comparison with v, v^ — 0, 

v =0, and v v = v. Then 
z 7 

kVSL . - L (tl + 2v + v 2 iff') 

p Dt c 2 \St 2 Sx St Sx 2 / 

With this relation vised, in equation (B3) the equation for the 
•velocity potential may be put in the form given in equation (4) 
of the analysis. 



and 
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APPENDIX C 


EVALUATION OF [ tt" 

\oy> 


y — *0 


In order to determine the limit of — ‘ as y— * 0, it is 

dy 

convenient to make use of the following substitution 


^ cos 6 + ^2 + ^1 


The expression for (ft (equa.tion (10) ) may "be written with the aid 
of the following relations (see equation ( 11 )): 


7 w w • « 

1 4 • L , 


j to sin e 
A r.« 


/*%. pTX 

0 * p Jq &(p 0, z +£ 0 cos 9)(f± + fg) ae d£ (C2) 


fl = f(t - T X ) = 


M(x - 5) L 8in 0 


ft - — — + 

V C P 2 


, \ f M(x - £) to sin 6 

t 2 , f(t - T ) , fit— i_2- - 

' ' V cp^ cp 
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By the rule for differentiation of a definite integral 

\ 



(C3) 


Make vine of the following relations : 

Sa to a (* + ip_ eos J) 


dy 

d^z + 

i 0 cos q) 

dy 



dA. 

dz dy 

COS 0 



a ( f i + f 2 ) 

5f l 

S(t - Tj) 

c)f g 

+ 3 

d(t " t 2 ) 


z K 1 ' + — 7 — 1 

Sy o(t - rV dy c(t - T g ) dy 


at 


-j > v 

- ■ “ ( f l ' f 2 ) — 8ln ® 

or. ‘ - cy 


Then, “by irrtegravion by parts , the next to the last integral in 

CiA / \ 

equation ( 03 ) 1 scones (with cos 0 d 6 - dv ; — (fn + fo) == u) 

cz oy ' / 


»St 


1a . x 

,l.; f + f \ sin 9 

Oz Air \ A ^ 


oy 


0 


*s 


fv. 


_ JL 

on 

d /„ 


1 fa 

i*-, + f^j — - - sin 9 


“ — 

--(f 

, - f 0 ] sin 0 cos 6 

Jo l> 

1 0 a z 2 

c(3 

dz 

dt ' " 

J- 


_ . .. d£> 

Jo Jc 

where the first term vanishes "because 3 in 9=0 at 9=0 and 0 
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, bL x ' 

Similarly, (with sin 0 d 0 = dv, — i - — - . — ~ 

■ oy 8t ' ' 

the last integral in equation (C3) "becomes 


- u) 




_A_ ^0 

L C P St 




cos 0 




+-y 


^1 r n 

Jo Jo 


1 SA -S V . 

— '■ — ^ — If', - f Q 1 sin & cos 0 
c6 Sz St V1 ' " 


A 5 2 / . p 

^?( f l + *2)' 008 e 
cp 2 St 2Vi ’ . 


d0 d£. 


where the first term vanishes "becaii.se ^1 = ^2 £ = 0 and 0 = « 

Then,"’ as y approaches zero from the positive sido;. there 
results in the. limit a contribution only from the first integral in 
equation (C3), 

= - l) A (x, 0 ; z) f(t) (C4) 

• : ^>1 ■" ‘ " ’ . • - 

Since — — changes sign as y changes sign, it follows that 

£y . , ' . • .• • • ■ 

as y approaches zero from the negative .side an equal and opposite 
result -is; obtained. ' 
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Figs. 1-3 



Figure 2.- Field of influence of spherical source moving at a constant 
supersonic velocity. Point (x, y, z) is affected at time t by two 
pulses originating at times t and t earlier. 
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Fig. 10 




